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ABSTRACT 

The formation of anionic species in the interstellar medium from interaction of linear 
molecules containing carbon, nitrogen and hydrogen as atomic components (polyynes) 
with free electrons in the environment is modelled via a quantum treatment of the 
collision dynamics. The ensuing integral cross sections are employed to obtain the 
corresponding attachment rates over a broad range of temperatures for the electrons. 
The calculations unequivocally show that a parametrization form often employed for 
such rates yields a broad range of values that turn out to be specific for each molecular 
species considered, thus excluding using a unique set for the whole class of polyynes. 
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1 INTRODUCTION 

In spite of the fact that the physical conditions in the 
Interstellar Medium (ISM) and in many of the plane- 
tary atmospheres are seemingly unfavourable to chemi- 
cal evolutions (e.g. low particle densities and fairly low 
temperatures) a very rich chemistry exists in such re- 
gions and has been follow e d by observations for a lon g 
time (|Snow fc McCalll l2006l : [Herbst fc van Dishoeckll2009l ). 
Among the recently discovered species, six different an- 
ionic molecules have been detected in the denser re- 
gions of the ISM, i.e. HC" (n= 4,6,8) and CnN" 
(n= 1, 3, 5) ( "McCarth y et al.l l2006l: 'Briinkcn etal] l2007l : 
ICernicharo et al.. .200711 ; lAgundez etal.: .2010. ) and further- 
more N-containing anions CnN~ (n= 1,3,5) w ere found in 
Titan's upper atmosphere l|Vuitton et al.ll2009l ). It is there- 
fore well established by now that several molecular anions 
have been found in the interstellar and planetary environ- 
ments, so that it is now crucial to understand as much as 
possible the ion chemistry of the species involved and the 
likely paths to the formation of such systems: the final rate 
constants for their formation are in fact important for the 
astrochemical modelling that tries to establish the most re- 
alistic reaction networks in that environment. 

Given the low densities of the mutual species present 
in the ISM, it is also a realistic possibility that such anionic 
molecules be formed in reaction with electrons, the latter be- 
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ing produced abundantly in the outer reg ions of the dense r 
clouds and in the planetary atmospheres ( Sakai et al.ll2007l '). 
Estimates of the electron temperatures in that environment 
suggest values that can go up to 1000 K or above. The reac- 
tions of such electrons with molecular species can therefore 
lead to the eventual formation of a negative ion through 
a variety of dynamical processes. For example, dissociative 
electron attachment (DEA) is a well known process 



M-H-f e" 



M" -hH, 



(1) 



whereby some of the bound atoms in the initial partner are 
detached after forming a stable residual anion, a process re- 
lated to the sign a nd the value of the electro n affinity (EA) of 
the M — H species (SebastiancUi et al. 201l'). Another option 
in the ISM environment is the radiative stabilization (RS) 
of the threshold resonances associated to the formation of 
a temporary negative ion (TNI) upon electron attachment 
close to zero energy 



M + e 



M' 



M" 



hv 



(2) 



a process which we have recently argued as possible for 
the formation of [C6H4]" l|Carelli et al.ll201ll ). The DEA 
processes are also known to be more likely to occur at 
low energies and to ch iefly involve a dehydrogenation path 
(|Graupner et al.ll2006l ) 



e" + HCCCN -+ CCCN" + H 



(3) 



which has, in fact, been suggested by our calculations as 
occu rring via resonant intermediates in several polyynes an- 
ions (|Sebastianelli fc Gianturcoll2012l ) . An additional feature 
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that should favor RS in competition with DEA is the oc- 
currence of virtual states and of zero-energy resona nces in 
unsaturated species like Ceo (jLucchese et al.l 1 19991 ) . Such 
special "compound states" are suggested to be very ef- 
ficient in "soaking up" t hreshold electron s, thereby lead- 
ing to anionic formation (|Field et al.ll200ll ') and have been 
found by our calculations to als o exist in anionic polyynes 
l|Sebastianelh fc Gianturcoll2012l ). 



2 CALCULATING ATTACHMENT RATES 

The question of having such variety of processes thus relates 
to scattering cross sections over the full range of low-energy 
regions where either DEA and RS processes are found to 
occur. They often reveal a broa d variety of resonances for 
metastable electron attachment jSebastianelli fc Gianturcol 
Such cross sections in turn lead to the corresponding 
rates given as 



A:TNi(r) = 



1/2 



(fcBT)2 



(4) 

In our present models it describes a convolution over electron 
temperatures of the integral, total cross sections ao{E) that 
include all metastable attachment processes for the given 
molecule over that range of energies, plus a sum over elas- 
tic and rotational excitation processes starting with \j) = 
states of the targets. In our procedure for the calculation 
of the relevan t cross sections - e.g. see for d etails our foun- 
dation paper: iLucchese fc Gianturcol l|l996l l - the channels 
associated with electronic excitations and vibrational exci- 
tations have been excluded to better model ISM conditions, 
although the entire quantum dynamics of the scattering of 
electrons from the polyynes has been car ri ed ou t using ab 
initio methods - see again: ICareUi et al.l (|201ll ). We also 
know, however, from previous analysis of polyatomic targets, 
that the inelastic cross sections contribute, at the relevant 
energy, a bout 10% of the integral c ross sections included in 
Eqn.Q - llrrera fc Gianturcol (|2005l ). 

One should also note that the ao computed cross sec- 
tions of Eqn.Q are obtained for fixed nuclear geometries 
of the target polyynes and therefore the needed nuclear dy- 
namical couplings with the electron motion, necessary to 
obtain DEA cross sections, are not included. However, they 
include direct and accurate evaluations from first principles 
of all the resonance formation states and all the virtual state 
effects. Thus, as we shall further explain below, the corre- 
sponding rates would be upper bound to the true rates since 
all integral cross sections are taken to lead the anionic sta- 
bilization. 

The whole electron attachment process, namely 



M(i?0+e" 



(5) 



can be divided into the primary attachment (labelled with 
/ctni) and a subsequent process (ka) that, via intramolec- 
ular vibrational rearrangements, can lead to a stable anion 
according to the following prototypical reaction: 



MiRfY 



(6) 



In the above equation ka refers either to the radiative stabi- 
lization (a — rs) or to the dissociative process {a = DEA), 



and 



[M-{Re)]* ^M{Rf)- +hi' 
[M-{R,)]* [M- A]- + A, 



both in competition with the auto-detachment process 



M{Re) + e- 



M" 



M{Rf) + e~ 



(7) 
(8) 

(9) 



In the above equations the changes in molecular initial ge- 
ometries, Ri, are indicated by the final geometries, Rf, or by 
the formation of an anion in its initial equilibrium geometry. 
Re- Note also that DEA forms a different molecular species 
[M- A]-. 

All the above processes are taken to occur via the for- 
mation of resonance states, transient negative ions (TNI) 
that play the role of gateway states to the final formation of 
all stable anions, as described in Eqn.(|SJ- 

Thus, the present calcula tions generate from first prin - 
ciple quantum-dynamics - see iLucchese fc Gianturcol (|l996l ) 
for all details - all the TNI states contributing to the rate 
of Eqn.(l5]), plus all the zero-energy resonances and virtual 
states that occur as the collision energy goes down to zero 
and which can further contribute to ^ea. 

One should also be reminded that the overall rate of 
electron attachment (EA), /cea, can be related to the rates 
of metastable anion formations, fcTNi, and to the autode- 
tachment and radiative stabilization rates, fcad and fcrs re- 
spectively (|Carellill20l3 ) 



fcEA 



^ad + ki 



(10) 



when anionic formation is assumed to be in a steady state 
(Hcrbst 1980). Since it is also assumed that fead S> fers one 
can further write 



^EA 



krs 
kad 



(11) 



Therefore the rates of electron attachment, ^ea, have a clear 
upper bound in the calculated fcxNi rates which use Eqn.Q 
unless fcad ~ fcrs. The present calculations for polyynes are 
thus providing from first principles electron attachment rate 
values which are, generally speaking, larger than the true 
rates of formation of anionic species. They are also shown 
here to exhibit a very marked dependence on the chemical 
features of the carbon-rich linear molecules examined in the 
present work. 



3 RESULTS AND DISCUSSION 

We have carried out detailed calculations of the integral 
cross sections, as in Eqn.@, over a broad range of colli- 
sion energies using the computational method well estab- 
lished in our group (jLucchese fc Gianturcd[l996l ') as already 
give n in detail in various pape r s on p olyyne specie s - e.g. 



io, 
; ISebastiancU i et al.l 



(|20ia ) 

It is also customary in evolutionary models of chemical 
networks in the ISM and in planetary atmospheres to rep- 
resent general two-body reactions via a parametric form of 



Electron- 



attachment rates for carbon-rich molecules 3 



Table 1. Fit coefficients for the polyynes presented in this paper, 
where a is in cm^ s~^, /3 is dimensionless, and 7 in K. Note that 



a{b): 



10''. See text and Egn. 11211 for further details. 
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/3 
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2.689(-10) 


0.486 


4.174 


NC2N 
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NC4N 


21.798(-10) 
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Figure 1. Computed fcTNl rates from Eqn.Q for three different 
carbon-rich hnear molecules without permanent dipoles. See text 
for further details. 



their c orresponding rate coefficients - e.g. see: IWoodall et al] 
|2003) 



( T 



V300K 



\ exp(-7/T) 



cm^s"^ , (12) 



where, in our case, T is the electron temperature Te- The 
present calculations can therefore be profitably represented 
by a set of three parameters: a, /3, 7, that can be associated 
to each polyyne species considered within the chemical net- 
work. 

The molecules which we consider correspond to primary 
molecules for which electron attachment rates can be calcu- 
lated: the symmetric polyynes (HC4H, NC2N, and NC4N) 
and the polar polyynes (HC3N and HC5N) as exemplary 
choices for the class of molecules which have been also ex- 
perimentally observed (see present introduction). 

The data reported by Fig[T]and Fig[2]show the temper- 
ature dependence of the fcTNi from the present calculations 
over a very broad range of Te'. from nearby zero to 5000 K. 

The data in both figures are then fitted using the cus- 
tomary expression of EQn. (|12p and the corresponding pa- 
rameter values are also reported (see also Tab[T]). 

The following features could be gleaned from these data: 

(i) There is a dramatic change in size when going from 
polar to non-polar molecular partners: the rates increase for 
the former by nearly two orders of magnitude over the whole 
range of T^; 

(ii) The polar targets also show a strong increase at low 
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Figure 2. Same as Fig[T]but for two polar polyynes. 



Te values, decreasing at larger temperatures; the opposite 
occurs for the symmetric polyynes, whose rates increase as 
Te increases; 

(iii) The parameters for the polar systems are fairly ro- 
bust, in the sense that change little with the number of car- 
bon atoms, a molecular index often used in astrochemical 
networks; 

(iv) The symmetric polyynes show instead a more marked 
dependence on the number of C atoms and also have larger 
rates with N atoms in the chains. 

The data insets of Fig[T] and Fig[2] also report the /ctni 
behaviour in the low-Te regimes, i.e. from threshold to 250 
K, in order to show more clearly the difference in behaviour 
among chemical partners. 

The marked increase of the rates near threshold and 
their drops at vanishing energies are again very clear for the 
polar molecules, as well as the larger sizes of their rates at 
low T. On the contrary, the symmetric polyynes show all 
a steady increase in size from threshold and overall sizes 
which are much smaller that those for their asymmetrical 
counterparts. 

In conclusion, the present quantum calculations of for- 
mation rates for electron attachment processes within a rep- 
resentative sampling of carbon-rich linear molecules show 
rather clearly that the parametrization of such two-body 
processes is markedly dependent on the chemical species be- 
ing considered. Although the length of the chains (i.e. the 
number of C atoms) plays some role, as does the presence 
of other heteroatoms like nitrogen, we found here that the 
strongest chemical discriminator is provided by the pres- 
ence of permanent electric dipole moments in the molecules 
in question. The asymmetric polyynes are seen here to yield 
much larger rates, especially for Te values from threshold to 
about 100 K, and to exhibit an increase at lower tempera- 
tures, the opposite behaviour from that shown by symmetric 
molecular partners. 

Such differences should therefore be taken into con- 
sideration when anionic formation rates are explicitly in- 
cluded in evolutionary chemical networks, where they can 
provide more detailed, chemically-controlled, final station- 
ary abundances for neutral and charged species. The latter 
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data would obviously be critical for a better c omparison with 
observations - e.g. see lAgundez et"aLl (|201(]| ). 



ACKNOWLEDGEMENTS 

The computational support from the CASPUR Consortium 
is gratefully acknowledged, as well as the financial support 
from the PRIN 2009 research network. One of us (T.G.) 
thanks the CINECA Consortium for the awarding of a post- 
doctoral grant during which this work was carried out. 



REFERENCES 

Agiindez, M., Cernicharo, J., Guelin, M., et al. 2010, A&A, 
517, L2 

Briinken, S., Gupta, H., Gottlieb, C. A., McCarthy, M. C., 
& Thaddeus, P. 2007, ApJL, 664, L43 

Carelli, F. 2012, "Molecular anions in circumstellar 
envelopes, interstellar clouds, and planetary atmo- 
spheres: quantum dynamics of formation and evolution", 
l^ttp : //arxiv . org/abs/1209 . 2827 

Carelli, F., Sebastianelli, F., Satta, M., & Gianturco, F. A. 
2011, Monthly Notices of the Royal Astronomical Society, 
415, 425 

Cernicharo, J., Guelin, M., Agiindez, M., et al. 2007, A&A, 
467, L37 

Field, D., Ziesel, J.-P., Lunt, S. L., et al. 2001, Journal of 
Physics B Atomic Molecular Physics, 34, 4371 

Graupner, K., Merrigan, T. L., Field, T. A., Youngs, 
T. G. A., & Marr, P. C. 2006, New Journal of Physics, 
8, 117 

Herbst, E. 1980, ApJ, 237, 462 

Herbst, E., & van Dishoeck, E. F. 2009, ARA&A, 47, 427 
Irrera, S., & Gianturco, F. A. 2005, New Journal of Physics, 
7, 1 

Lucchese, R. R., Gianturco, F., & Sauna, N. 1999, Chemical 
Physics Letters, 305, 413 

Lucchese, R. R., & Gianturco, F. A. 1996, International 
Reviews in Physical Chemistry, 15, 429 

McCarthy, M. C, Gottlieb, C. A., Gupta, H., & Thaddeus, 
P. 2006, ApJL, 652, L141 

Sakai, N., Sakai, T., Osamura, Y., & Yamamoto, S. 2007, 
ApJL, 667, L65 

Sebastianelh, F., CareUi, F., & Gianturco, F. 2012, 
Chemical Physics, 398, 199 , chemical Physics of Low- 
Temperature Plasmas (in honour of Prof Mario Capitelli) 

Sebastianelh, F., CareUi, F., & Gianturco, F. A. 2011, The 
Journal of Physical Chemistry A, 115, 11531 

Sebastianelli, F., & Gianturco, F. A. 2012, European Phys- 
ical Journal D, 66, 41 

Snow, T. P., & McCaU, B. J. 2006, ARA&A, 44, 367 

Vuitton, v., Lavvas, P., Yelle, R. V., et al. 2009, 
Planet. Space Sci., 57, 1558 

Woodall, J., Agiindez, M., Markwick-Kemper, A. J., & Mil- 
lar, T. J. 2007, A&A, 466, 1197 

This paper has been typeset from a IJjX/ I^TpTjX file prepared 
by the author. 



